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Abstract 

We demonstrate that energy can be transferred very rapidly from the bulk to 
the brane through parametric resonance. In a simple realization, we consider a 
massless bulk field that interacts only with fields localized on the brane. We assume 
an initial field configuration that has the form of a wave-packet moving towards the 
brane. During the reflection of the wave-packet by the brane the localized fields can 
be excited through parametric resonance. The mechanism is also applicable to bulk 
fields with a potential. The rapid energy transfer can have important cosmological 
and astrophysical implications. 



PACS numbers: ll.10.Kk, 04.50, 98.80.Cq 



Introduction: The idea that our Universe is a defect embedded in a higher-dimensional 
bulk space is old [1]. The Randall-Sundrum model [2] provides a realization of this idea. 
The defect is a four-dimensional hypersurface (a 3-brane) in a five-dimensional bulk with 
negative cosmological constant (AdS space). The geometry is non-trivial (warped) along 
the fourth spatial dimension, so that an effective localization of low-energy gravity takes 
place near the brane. The matter is assumed to be concentrated only on the brane. 
For low matter densities the cosmological evolution on the brane becomes typical of a 
Friedmann- Robertson- Walker Universe [3, 4]. 

Of particular interest in the scenario of a brane Universe is the possibility of exchange 
of energy between the defect and the bulk space. For this to occur the bulk must contain 
some matter component in addition to the negative cosmological constant. The presence 
of matter in the bulk affects the expansion on the brane. For an observer on the brane the 
modification can be attributed to "mirage" matter components [5]. The energy exchange 
can modify the cosmological evolution as well [6] . In particular, the energy transfer from 
the bulk to the brane can lead to periods of accelerated expansion on the brane [7]. 

In this letter we concentrate on the mechanism that can lead to a rapid energy transfer 
from the bulk to the brane. We omit the gravitational effects and study a simplified 
system that consists of a flat five-dimensional bulk space with a four-dimensional brane 
as its boundary. In the bulk we assume the presence of a field that interacts only with 
fields localized on the brane. Our aim is to demonstrate that energy can be transferred 
very efficiently from the bulk to the brane fields through a process that resembles very 
strongly the parametric resonance [8, 9]. The incorporation of gravitational effects is 
technically complicated and will be left for a future publication. 



The model: We consider a massless bulk field <p that interacts with a field x localized 
on the brane. We neglect gravitational effects and assume that the system is described 
by the action 

S = J d A xdy[d M ^d M ^ + 5{y){d» X d»X-V{x^))\, (1) 

where M — 0, 1, 2, 3, 4 and // = 0, 1, 2, 3. In order to keep our discussion simple we identify 
the half-space y < with the half-space y > 0, in complete analogy to ref. [2]. In this way 
the brane forms the boundary of the bulk space. We assume that the potential V(x, <P) 
has the form 

V{xA) = \rn\ 2 + \gx 2 <t> 2 - (2) 
We expand the quantum field x as 

X(t, *) = / + *lr k (t)e fk£ ) , (3) 

where the mode functions fk satisfy fkft~ fkfk = h and the number distribution function 
is n k = (aj-Ofc). The equation of motion of x gives 

f k (t) + (k 2 + m 2 + g<j>\t, 0)) f k {t) = 0, (4) 



2 



where a dot denotes a time derivative. We also have 

/d k* 
J^n(t)h(t)(l + 2n k ). (5) 

We treat the field <fi classically. We also make the simplifying assumption that the 
configuration of depends only on the time t and the fourth spatial coordinate y. For 
this assumption to be consistent with the 0-field evolution we replace x 2 (^ ^0 with (x 2 )(^) 
in the equation of motion of 0. This is the Hartree approximation [10] that is usually 
employed in the study of parametric resonance [9] . The 0-field equation of motion becomes 

0(t, y) - <j>"(t, y) + g 5{y) 0(y, t) ( X 2 ) (t) = 0, (6) 

where a prime denotes a y-derivative. Integrating this equation with respect to y in the 
interval [— e, e] and taking e — > 0, we find 

0'(t,O+) = -<f>'(t,0-) = ±g<f>(t,0) ( X 2 )(t). (7) 
The solution of eq. (6) is equivalent to the solution of the bulk equation 

4>(t,y)-<f>"(t,y) = (8) 

with the boundary condition (7). 

In the following we solve numerically the system of equations (4), (5), (7), (8) in the 
region y > 0. The solution for negative y is <f>(t,—y) = <f)(t,y). The initial configuration 
0(0, y) is arbitrary. We choose it to have the form of a wave-packet that moves towards 
the brane. For the mode functions of the x-field for t — > 0, when 0(£,O) ~ 0, we assume 
the vacuum expressions f k (t) = e~ lWkt / \J2ojk, with uj k = \Jk 2 + m 2 . This gives fk{0) = 
1/ \/2ui kl fk(0) = —ioj k / \J2uj k . The evolution of f k {t) can be attributed to the creation of 
X-particles on the brane through the time-varying 0-field. The time-dependent number 
distribution function can be calculated through a Bogoliubov transformation [10]. It is 
approximately 
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fk(t) 
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n k (t) = ^u k (t) \\f k (t)\ 2 + J J^J (2n k + 1) - ^ (9) 

with u k {t) = \jk 2 + m 2 + g<fi 2 (t, 0). The above expression has corrections ~ u) k that we 
have neglected by assuming that io k {t) w k (t)- In the limit t — > oo, after the 0-wave- 
packet has been reflected by the brane and <f>(t, 0) ~ 0, our approximation becomes exact. 
This means that the final number of x-particles is computed exactly, even though the 
expression (9) is approximate when <f>(t, 0) ^ 0. The initial number distribution function 
n k can be chosen according to our assumptions for the initial state of the physical system. 
For example, it can have the form of a thermal distribution [10]. We make the simpler 
assumption that initially there is no significant number density of x-particles, so that 
n k = 0. 
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Several approximations that we have made in deriving the evolution equations for 
the fields are not expected to affect the qualitative physical behaviour. For example, a 
term ~ x 4 m the potential of eq. (2) has been taken into account in previous studies of 
parametric resonance, without altering significantly the essential behaviour. The Hartree 
approximation becomes exact in the large- N limit, if the field \ is assumed to have N 
components and the Lagrangian to be invariant under an O(N) symmetry. Finally, the 
possible presence of a potential for the bulk field is not essential either. As we shall see 
in the following, the necessary requirement is that the bulk field must be oscillatory at 
the location of the brane. Any wave solution that is reflected by the brane is capable of 
transferring energy onto it through parametric resonance. 



An analytical solution: In order to have some intuitive understanding of the solutions 
of interest we derive in this section an analytical solution of eqs. (4), (5), (7), (8). We 
make the assumption 0(t, 0) = (p cosat. Then, eq. (4) can be put in the form of the 
Mathieu equation 

' 2 Tk 



d 2 fk 



+ (A k -2qcos2z)f k = 0, (10) 



dz 2 

with z = at — 7r/2, A k = (k 2 + m 2 + g<j)l/2)/a 2 , q = g<fil/2a 2 . An infinite number of 
resonance bands exist in k, within which the amplitude of f k grows exponentially fast 
with time. As a result, (x 2 )(t) grows fast. 

The form of <j)(t, y) that is consistent with the above solution for x can be determined 
analytically. It is given by the solution of the wave equation (8) for y > 0, with the 
"initial" conditions given for y = 0. Essentially the role of t and y is reversed compared 
to the standard case. The most general solution of eq. (8) is 



rt+y 

a(t + y) + a(t-y)+ / b(t')dt' 

Jt-y 



(11) 



where the functions a(t), b(t) are determined through the boundary conditions at y — 0. 
In particular, a(t) = 4>(t, 0) = 0o cos at according to our assumption. The function b(t) is 
defined through eq. (7), that gives 6(f) = <j)'(t,0) = gcf)(t,0)(x 2 }(t)/2. 

For a given time t the field (f>(t, y) has an oscillatory form with an envelope that grows 
as a function of y. As time increases the slope of the envelope grows. This behaviour 
is a consequence of our assumption that has a constant amplitude 0o at y = 0. The 
growth of (x 2 ){t) induces a backreaction on the evolution of <fi that tends to reduce its 
amplitude. This effect can be compensated only if the amplitude of the bulk oscillations 
of increases with time. 



The numerical solution: The analytical solution indicates that the rapid transfer of 
energy from the bulk to the brane through parametric resonance is possible. However, 
the solution was constructed by assuming the form of the 0-field at the location of the 
brane and deducing its form in the bulk. Here we would like to start by defining initial 
conditions for the system, such that the deviation of from zero is significant only away 
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Figure 1: The form of the bulk field at various times. 



from the brane. Subsequently, the 0-field perturbation moves close to the brane and leads 
to the excitation of the x field. This calculation cannot be performed analytically because 
of the difficulty in implementing the boundary condition (7). For this reason we solve the 
system of equations (4), (5), (7), (8) numerically. 

For early times we consider a configuration for the 0-field given by 

<P{t, y) = 0o exp — ^° 2 + ^ ^ cos [a(y - y + t)} . (12) 

It corresponds to a wave-packet of fundamental frequency a and width a, centered at 
y = i/o — t, that moves towards the brane located at y — with the speed of light. In the 
following we study the evolution of a wave-packet with O = 10, y — 5, a — 1, a — 20. 
We use the values m 2 — 1, g — 10 for the parameters in the potential of eq.(2). Essentially, 
we express all dimensionful quantities in units of m. We also assume a negligible initial 
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Figure 2: The real part of the mode function fk of the brane field at various times. 



particle number density: n& = 0. The total number density of x-particles and their energy 
density are given by 

N(t) f d 3 k „, E(t) f d 3 k , , 

~V = J TZnfi ®' T = /(2i)i"> (i »'"<* (13) 

with rik(t) given by eq. (9). The momentum integrations in these equations, as well as 
in eq. (5), have ultraviolet divergences. We regulate them through an explicit cutoff 
A = 10. Moreover, the time-evolution of the x-field is expected to display the well-known 
ultraviolet problem of classical systems. (The Hartree approximation in the leading order 
leads to essentially classical evolution.) The high-frequency modes become excited at 
sufficiently large times and the ultraviolet contributions dominate in the expressions for 
the various observables. In our simplified study we omit all x-field self-interaction terms. 
As a result, the ultraviolet modes can be excited only during the finite interval that 
the wave-packet interacts with the brane. For our choice of parameters, no significant 
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Figure 3: The number distribution function on the brane at various times. 



excitation of these modes takes place. 

In fig. I we display the evolution of the wave-packet. At t — it is located away from 
the brane. Subsequently it moves towards the brane and interacts with it during the time 
interval 3 < t < 7. Eventually it gets reflected with a distorted profile. At times t ~ 5 the 
maximum of <p(t, 0) is ~ 20o- The reason is that, because of our assumption of reflection 
symmetry around y — 0, there is another wave-packet interacting with the brane from 
the left. 

In fig. 2 we depict the evolution of the mode functions fk(t) of the %-field. We observe 
that they increase strongly in the region of small k through parametric resonance. The 
mode functions with k near the cutoff A = 10 are not excited significantly during the time 
interval that we consider. This implies that we can trust the predicted physical behaviour. 

In fig. 3 we display the number distribution function rik(t) as given by eq. (9). We 
observe that during the time interval 3 < t < 7 the low-frequency modes are excited. The 



7 




Figure 4: The total number density of x-particles, the energy density carried by them and 
the dispersion of the x-field, as a function of time. 



quantity nk(t) displays some oscillatory behaviour with time, but on the average it grows 
fast. For t > 7 the wave-packet moves away from the brane and the distribution function 
becomes constant. This is a consequence of the omission of any x-field self- interactions 
in our simplified treatment. 

In fig. 4 we plot the total number density of x-particles, the energy density carried by 
them and the dispersion of the x-field. All these quantities display a rapid increase in the 
time interval 3 < t < 7, while they remain constant at later times when the wave-packet 
has moved away from the brane. 

The intensity of the parametric resonance depends on the parameters of the model. We 
can use the Mathieu equation (10) in order to obtain some intuition. The amplification 
is more pronounced if the system is in the region of broad parametric resonance. This 
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requires gcfil/m 2 ^ 1, so that A k <: q. The parameters we chose satisfy this constraint. 
The time dependence of n k (t) in fig. 3 is typical of a system in the region of broad 
parametric resonance [9]. 

The resonance effects could be further enhanced if the effective mass became zero for a 
certain value of \- This would require m 2 < and would imply that the 0-field determines 
the x expectation value. We have used a simpler model, in which the bulk plays only a 
minor role in the determination of the vacuum state of the brane fields, while it can affect 
their dynamical evolution. 

Despite the similarity we described above, a careful analysis indicates that the Mathieu 
equation is not strictly relevant if the bulk field has the form of a wave-packet. The 
reason is that the function <p(t, 0) that appears in eq. (4) describes an oscillation with an 
amplitude that is non-zero only within a finite time interval. The Fourier transform of this 
function is non-zero over a large range of frequencies, and not for a unique frequency as 
in the Mathieu equation. We expect a qualitative behaviour similar to the one predicted 
by the Lame equation. In agreement with this conclusion, we do not find any parameter 
range that resembles the region of narrow parametric resonance. The numerical analysis 
of the range g(pl/m 2 ^ 1 leads to an evolution very similar to the one depicted in figs. 1-4. 
The only difference is that the number density of the produced x-particles and the total 
energy trasferred to the brane fall very fast with decreasing gcf^/m 2 . When gcfl/m 2 <C 1 
the resonance effects become negligible. Similarly, significant excitation of the x _m °des 
is possible only if the fundamental frequency a of the <fi wave-packet is not much smaller 
than the mass m. 

Cosmological and astrophysical implications: The rapid transfer of energy from 
the bulk to the brane can affect the cosmological evolution. If the x-field is coupled to 
lighter brane particles the energy can be transferred to them through the decay of the 
X-modes. This would result in a period like preheating [8, 9]. The two mechanisms are 
so similar that the process we have described could be characterized as preheating by a 
bulk source. There are several advantageous aspects for this type of preheating: 

a) The form of the source is less constrained than in the conventional scenario. In partic- 
ular, the field <fi need not be homogeneous. In the example we discussed it has the form 
of a wave-packet that depends only on the extra spatial coordinate. There is no need for 
a specially chosen potential for 0. In our example the potential was assumed to be zero. 

b) If is not the field responsible for the brane inflation, the time of preheating and the 
amount of the released energy are not constrained. In our example, for gcfl/m 2 ^ 1 the 
energy density in the x~ m odes exceeds m 4 . For large m significant amounts of energy can 
be deposited on the brane. 

On the other hand, this scenario requires an additional non-compact spatial dimension. 
A picture consistent with observations can be obtained only by generalizing the Randall- 
Sundrum model [2] in a cosmological context. The graviton localization near the brane 
leads to a low-energy cosmological evolution typical of a four-dimensional Universe [3, 4] . 
The effect of the energy influx on the cosmological evolution can be more complicated than 
the mere increase of the brane energy density. It is possible that an era of cosmological 
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acceleration can take place [7]. A detailed investigation of this issue must include the 
solution of the Einstein equations for the generalized warped geometry of the Randall- 
Sundrum model [4]. A numerical study of such a system is technically difficult. The 
only existing study has considered a two-brane scenario without brane fields [11]. The 
numerical investigation of the scenario with brane fields is in progress. 

Another manifestation of the mechanism we described could be related to the energy 
spectrum of high energy cosmic rays. It has been argued that these could be produced 
in the decays of long-lived heavy particles of cosmological origin [12, 13]. Our mechanism 
provides an alternative possibility: The decaying heavy particles could correspond to the 
X-modes that are produced during the energy transfer from the bulk to the brane. In order 
to provide an explanation for the highest energy cosmic rays, the mass of the x-particles 
must be m ^ 10 13 GeV [12, 13]. The production of so heavy particles is possible only if 
the fundamental frequency a of the wave-packet is of the same order of magnitude. 

The total amount of transferred energy is controlled by the combination gfo, where g 
is the strength of the coupling between the bulk and brane fields and 4>o the maximum 
amplitude of the wave-packet. The total amount of energy stored in x-particles drops 
very fast for g<p <C m, with m the mass of the x-field. In this range the resonance effects 
become negligible. Moreover, the dilution of the energy density through the possible 
accelerated expansion on the brane because of the energy influx has not been taken into 
account in our discussion. The combination of the two factors indicates that the total 
amount of transferred energy can be much smaller than m 4 , so as not to disturb the 
conventional cosmological picture, even though the energy will be distributed in heavy 
particles of mass m. The quantitative analysis of these issues will be possible only through 
the numerical solution of the Einstein equations. 

As a final comment we mention that the role of the bulk field can be played by the 
gravitational sector of the theory. It is possible that a perturbation of the bulk metric 
that has the form of a gravitational wave interacting with the brane can result in the 
excitation of heavy brane particles. This is an effect opposite to the production of bulk 
gravitons by the brane matter [14]. 
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